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Abstract 

TmGa3 (AuCu3 structure) undergoes two phase transitions, an antiferroquadrupolar transition 
at ~ 4.29 K and long-range antiferromagnetic ordering at ~ 4.26 K. Due to the close vicinity of 
the two phase transitions, TmGa3 offers an interesting system to study the interplay of charge and 
magnetic degrees of freedom. In order to understand this interplay we have performed inelastic 
neutron scattering experiments on TmGa3 in the paramagnetic regime (T > 5 K) to redetermine 
the crystal electric field level scheme. By fitting our spectra at various temperatures we obtain a 
new crystal field level scheme with Lea, Leask and Wolf parameters of xllw = -0.44(2) and W 
= -0.222(2) K. The total crystal field splitting at 5K amounts to ~ 2.3 meV, about an order of 
magnitude less than found previously, but in good agreement with the splitting extrapolated from 
the related ErGa3 system. 

Our analysis yields a I^ singlet as the crystal field ground state followed closely by a (nonmag- 

(2) 

netic) Ti singlet at 0.009 meV. The next excited states are a I5 triplet at ~0.5 meV, which is 
almost degenerate to a I^ doublet. This level scheme is adverse to previous findings. Subsequent 
analysis of the magnetisation along several crystallographic directions and the temperature depen- 
dant susceptibility as well as of the magnetic contribution to the specific heat are consistent with 
our new crystal field parameters. Implications for the antiferroquadrupolar and the antiferromag- 
netic transition are discussed. 
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I. INTRODUCTION 



In 4/ electron systems the determination of the crystal electric field (CEF) level scheme 
is an important first step to understanding the rich variety of physical phenomena that 
rare earth compounds exhibit. In 4/ systems the coupling between the spin and orbital 
angular momentum is stronger than the crystal field due to the fact that the 4/ orbital 
lies deep within the ion core and the other occupied orbitals screen out the potential of 
the surrounding ions. To a very good approximation the magnetic properties of the free 
ion can therefore be characterized by the total angular momentum J. CEF interactions, in 
general, lead to a splitting of the 2J+1 manifold into CEF states. Their degeneracy and 
energy splitting is determined by the symmetry of the electric charge distribution and the 
magnitude of the interaction. 

Recently, there has been a renewed interest in systems which undergo antiferroquadrupo- 
lar (AFQ) phase transitions originating from the interaction of low lying quadrupole active 
CEF levels. CeBg was the first compound to be unambiguously defined as having an AFQ 
character by INS experiments. 1 The phase diagram of CeB 6 shows a huge enhancement of the 
AFQ phase with increase of magnetic field aligned along [001]. More recently, a modulated 
AFQ structure was identified in the intermetallic compound PrPb3 by INS experiments- A 
non magnetic Ti ground state was found in the heavy fermion superconductor PrOs4Sbi2 
which also undergoes a field induced O yz type AFQ ordering.-^ The heavy fermion sys- 
tem PrFe4Pi2 shows a large increase of the lattice thermal conductivity at the quadrupolar 
transition.- In UPd3 a Q zx quadrupolar order parameter was determined to describe the 
AFQ phased 

Quite a number of systems with AFQ transitions were found among Tm based com- 
pounds. TmZn (Refi7), TmCu (RefQ) and TmGa3 (Refjy) exhibit an AFQ order in the 
paramagnetic phase, followed closely by an antiferromagnetic (AFM) transition. TmCd 
(Refj7|), on the other hand, shows only a transition of AFQ nature, with no experimentally 
verified AFM transition down to 30 mK. TmTe, a magnetic semiconductor which contains di- 
valent Tm 2+ ions, has been studied by INS, in which a field induced magnetic superstructure 
was seen due to AFQ ordering at 1.8 K. 10 TmGa3 was intensively investigated due to its com- 
plex low temperature phase diagram. The phase diagram of TmGa 3 was first investigated 
by Czopnik et al. by specific heat, thermal expansion and susceptibility measurements.- 
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They found two close lying transitions at 4.26 K and 4.29 K. By a comparison with TmZn 
and TmCd, the first transition at 4.29 K was attributed to a structural transition driven by 
quadrupolar pair interactions, followed by an antiferromagnetic transition at 4.26 K. This 
interesting interplay of charge and magnetic degrees of freedom refocused our attention on 
TmGa3. TmGa3 crystallizes with the AuCu3 type cubic structure, space group Pm3m. The 
trivalent state, Tm 3+ , with electronic configuration 4/ 12 gives a total angular momentum J 
= 6. A first INS study carried out by Morin et al.— concluded the crystal field parameters 
in the Lea, Leask and Wolf (LLW) scheme (Ref 
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to be x = -0.32 and W = 1.03 K. These 



give a triplet as the CEF ground state separated by ~2.5 meV from a 1^ state as the 
first excited state. A Ig triplet CEF ground state was reported to be consistent with the 
temperature dependance of the magnetic susceptibility. A 169 Tm Mossbauer spectroscopy 
study was used to determine the quadrupolar term in the CEF Hamiltonian.— By fixing W 
= 1.0 K, 2- llw = -0.38 and using an iterative procedure the quadrupolar CEF term B 2 = 
-0.30 K was obtained. 

We decided to extend these preceding experiments and to search for possible splittings 
of the ground state, as well as follow the temperature dependance of the low crystal 
field levels near the AFQ transition to see whether there are any signatures of the AFQ 
transition in the CEF spectra. However, the results of our inelastic crystal field spectroscopy 
experiments reveal a much different CEF level scheme and overall splitting of the J=6 
manifold in TmGa3 that imply a reinterpretation of the low temperature behavior of TmGa3. 



II. THEORETICAL BACKGROUND 



A prerequisite in the analysis of the magnetic properties of TmGa3 is the determination 
of the CEF ground state of the Tm 3+ ions and the energy separation to the excited levels. 
For a cubic system the Hamiltonian of the 2 J+l manifold of the total angular momentum J 
due to a crystal field interaction of cubic symmetry is conveniently written in the operator 
equivalent forrn^ 



Hcef = B 4 (0° A + 50 4 4 ) + £? 6 (0° - 210(1) 



(1) 



where are the Stevens operators tabulated e.g. in Ref. 



14J . The coefficients and 



Bq are parameters that measure the respective components in the multipole expansion of 



4 



the CEF potential. With the relations 



B A = A A <r A >< J\\(3\\J> (2) 

and 

B 6 = A 6 < r 6 >< J||7||J >, (3) 

B A and B§ can be further decomposed into a product of the coefficients A A and 
sometimes called geometrical factors^, the mean fourth and sixth powers of the radial part 
of the wave functions of the 4/ electrons, < r 4 > and < r 6 > and the factors < J\\/3 , 7|| J > 
which are listed e.g. in Ref. 141 ]. 

In order to conveniently tabulate the normalized eigenvectors and eigenvalues for a given J 
manifold, Lea, Leask and Wolf suggested a transformation of B A and Bq into the parameters 
^llw and W by using the following relations^ 



5 4 F(4) = Wx LLW 
B 6 F(Q) = W(l-\x L1M \) 



(4) 
(5) 

where F(4) and F(6) are numbers given in Ref. [12|. Via this transformation the whole 
range covered by B A and B$ is mapped onto the the variable xllw limited to the interval -1 
< x < +1 and the energy scale factor W. With this transformation eq.([T]) now becomes 



H, 



CEF 



w 



(^iy) + (1 - |XLLwl) (^);j 



(6) 



with 



o A = ol + w\ 

6 = 6 ° - 210 6 4 



(7) 
(8) 



Empirically it was found that for a given system of compounds with the same crystal 



15J6 



structure the factors A A and A 6 vary only a little across the rare eart 
sequently, by using tabulated values for < r 4 > and < r 6 > (RefsJl? 
< J||/3||J > and 7 j = < J||7| J > one can estimate unknown parameters B A and Bq from 



ns series.— 1 ^ Con- 
18l ) and for (3j = 



known CEF splitting found for another rare earth ion in the series of isotypic compounds. 



Morin et al., from the analysis of their INS spectra (see above), found LLW parameters 
x LLW =-0.32 ± 0.02 and W=1.03 ± 0.03 K for TmGa 3 , corresponding to A A < r 4 > = -34 ± 
4 K and A 6 <r 6 > = -17±2K with a Ig triplet as the CEF ground state and an overall 
splitting to the uppermost T 2 singlet of 17.5 meV.— 

A detailed analysis of the crystal field splitting of the Er 3+ manifold in the system ErGa 3 
(isotypic to TmGas) was carried out by Murasik et al. by inelastic neutron spectroscopy 
and by the magnetic contributions to the heat capacity and the magnetization 19 . These 
experiments resulted in LLW parameters for Er 3+ of £llw — +0.195 and W = +0.022 meV 
or equivalently, using eqs. (jSJ) and (jSJ), gives A 4 (Er 3+ ) = +12.3 K and A 6 (Er 3+ ) = +3.85 K. 
Using the numerical values of < J\\/3\\ J > and < J||7|| J > for Er 3+ and Tm 3+ and also the 
values for < r 4 > and < r 6 > of Er 3+ and Tm 3+ respectively^ 1 ^, we find a large discrepancy 
in the magnitude and sign, for the coefficients A4 and Aq for Er 3+ in ErGa3 to those of Tm 3+ 
in TmGa3, as derived from the result of Morin et al.. 

TmGa3 undergoes quadrupolar and antiferromagnetic ordering below ~4.3 K and ~4.2 
K, respectively (see above). Due to magnetoelastic coupling, associated to these transitions 
is a distortion of the lattice which leads to a symmetry lowering from cubic symmetry.- 1 ^ 
This distortion induces a further splitting of the CEF states that can be described by an 
additional term added to the Hamiltonian in eq. © 



Quadrupolar-quadrupolar interaction between a pair of rare earth ions labeled i and j is 
usually described by the Hamiltoniani 4 -^^ 



QP — 



BlO% + B\0\. 



(9) 



H QQ = A [AOl Ol - 16 [OH Q£ + 2 l 0+ 1 ) + (0+ 2 2 f + 2 f 0+ 2 )] 



(10) 



where the operators O k 



are defined in the usual way by 



0° 2 = 3J 2 Z - J(J+1), 



(11) 




(J Z J± + J± J z ) , 



(12) 



and 



,±2 




(13) 



6 



The coefficient A is given by 



3e 2 < rf >< r\ >< Jj||a||Jj >< J,-||a||Jj > 

A = ; (14) 

with < rf > being the mean square radius of the 4/ ions on each site and < J||a||J > 



the coefficient listed e.g. in Ref. 
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R is the distance between the ions and is an 
effective dielectric constant for the pair of ions.— Knowing < r 2 > and e^, the coefficient 
A can, in principle, be calculated. However, is very sensitive to shielding effects due to 
the intervening ions and the conduction electrons in the case of metals. In general, A must 
therefore be treated as an unknown parameter to be obtained from experiment. 

First order quadrupolar effects for Tm 3+ may occur for the non-Kramers doublet T 3 with 
the two states having a quadrupole moment of opposite sign but the same magnitude and 
for the triplet states T 4 and T 5 . The two singlets I\ and T 2 show only second and higher 
order quadrupolar effects. 



III. EXPERIMENTAL 



All samples were grown from a Ga flux in quartz ampoules with a procedure described 
elsewhere^ and characterized by X-ray diffraction to ensure phase purity. They were found 
to have the correct AuCu3 structure with no impurity reflections above the noise level. The 
composition of single crystals was checked by electron microprobe analysis and found to 
have the correct atomic ratios. 

Inelastic neutron CEF spectroscopy was performed on the neutron time-of-flight spec- 
trometers IN4, IN6 and the triple axis spectrometer INI 2 at the Institut Laue-Langevin 
(Grenoble). IN4 operates in the thermal neutron energy range 10-100 meV. A suitable en- 
ergy is selected from the thermal spectrum with a crystal monochromator. The cold neutron 
time-of-flight spectrometer IN6 provides quasielastic and inelastic spectra for incident wave- 
lengths in the range of 4 to 6 A. A graphite monochromator delivers the four wavelengths 
4.1, 4.6, 5.1, and 5.9 A. The second order reflection from the graphite monochromator is 
removed by a beryllium-filter cooled to liquid nitrogen temperature. The elastic energy 
resolution at 4.1 A is 170 /xeV and at 5.9 A, 50 /xeV. Inelastic neutron scattering measure- 
ments on a polycrystalline sample of Tmo.1Luo.9Ga3 were carried out on the cold neutron 
triple-axis spectrometer IN12. The analysis of the INS was carried out using the program 
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Mcphasa^ to calculate at a given temperature the intensities and energies of the allowed 
CEF transitions for a set of LLW parameters xllw and W. 

In addition, we have performed specific heat measurements on single crystals and poly- 
crystalline samples of TmGa3 and polycrystalline samples of Lui-^Tm^Gaa (x=0.05, 0.01) 
and on the isostructural nonmagnetic compound LuGa 3 in a PPMS (Quantum Design) 
relaxation-type calorimeter. The sample pieces were fixed with a minute amount of Apiezon 
N grease on a sapphire platform which also carried the heater and a calibrated temperature 
sensor. The addenda heat capacities of the calorimeter and the grease were determined in 
separate runs and subtracted. Magnetic susceptibility measurements were performed us- 
ing a MPMS magnetometer (Quantum Design) in the temperature 1.9K < T < 350K and 
0<T< B < 7T. 

Electron paramagnetic resonance (EPR) spectroscopy on Er 3+ in polycrystalline samples 
with composition Lui-^Er^Gas (x=0.005, 0.01, 0.03) was performed using a Bruker X-band 
EPR spectrometer [y ~9.3 GHz). The derivative of the resonance absorption dP a bs/dH 
is obtained using field modulation (u=100 kHz) with standard lock-in detection technique. 
Coarse powdered samples were filled in quartz ampoules and fixed in paraffin for mea- 
surements below room temperature, which was provided by a continuous-flow He cryostat 
(Oxford instruments). 

IV. RESULTS AND DISCUSSION 
A. EPR on Lui x Er x Ga 3 

In order to test whether the finding of a T 7 doublet CEF ground state for the J = 15/2 
manifold of Er 3+ in ErGa3 can be confirmed, we carried out an EPR study on Lui-^Er^Gas 
(a;=0.005, 0.01, 0.03) in the temperature range 4 K < T 

lesssim25 K. A T 7 Kramers doublet is expected to exhibit an isotropic EPR line at a 
resonance field corresponding to a g factor of g{T 7 ) ~ 6.8. The resonance field of the T 7 
doublet can be well distinguished from that of the r 6 doublet which has a g factor of g(T 6 ) « 

(1—3) 

6.0. EPR resonances of the three Tg quadruplets resulting from the CEF splitting of the J 
= 15/2 manifold are angular dependent and in a polycrystalline sample will generally result 
in very broad smeared resonance lines. Excited T g CEF states at an energy A ~ ksT open 



S 



additional relaxation channels via Hirst and Orbach processes which give an exponential 
increase of the linewidth, oc exp(—A/kBT)) in addition to the normal Korringa relaxation, 

Murasik et al.'s CEF scheme provides a Ig quadruplet as the first excited state at an 
energy of ~2.7 meV.~ Measurements of the temperature dependence of the EPR linewidth 
above liquid He temperature should be able to reveal relaxation via this state and provide 
an additional support for the validity of the proposed CEF level scheme of ErGa3. 

As a characteristic spectrum, Fig. [1] displays the EPR resonance of Lu .995Er o5Ga 3 at 
a frequency of ~9.3 GHz. At low temperatures we observe an asymmetric resonance line 
at resonance fields of ~ 978(1) Oe with typical linewidths of ~100 Oe. With increasing 
temperature the resonance rapidly broadens and it is not detectable any more above ~25 
K. As is typical for EPR of localized moments in metals and intermetallic compounds, the 
resonance is asymmetric due to the mixtures of an absorption and dispersion signal due to 
the skin effect. Characteristic shoulders, especially resolved for small Er concentrations in 
the low field wing of the central line are due to hyperfine satellites from the isotope 167 Er (/ 
= 7/2). Their resonance fields H res (m = —7/2, ...,7/2) shift with respect to the resonance 
field of the central line, H , observed for all other isotopes with 1=0 according to: 26 

167 /t 167/12 

H res {m) = H - 167 Am[l + (— -) 2 ] - —-[1(1 + 1) - m 2 ] (15) 
where 167 A is the hyperfine constant for 167 Er. 

A superposition of the central line and the hyperfine satellites, all with the same linewidth, 
and an absorption/dispersion ratio (characterizing the linewidth asymmetry) of typically ~1 
was found to fit the spectra very well (cf. Fig. [Q). 

The resonance fields of the central lines is independent of the temperature within error 
bars and corresponds to a g factor of 

g = 6.798(4). (16) 

The g factor and hyperfine constant 167 A are somewhat larger than the respective quan- 
tities observed for the Er EPR of a T-j CEF doublet in insulators (6.75 - 6.76 and 73 - 74 
G).— These shifts can be understood as due to the exchange interaction of the localized Er 
moment with the conduction electrons.— 1 ^ 
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FIG. 1: (color online) o EPR resonance line of Er 3+ in the sample Luo.ggsEro.oosGas at 4 K at a 
microwave frequency of 9.3034 GHz. The (red) solid line is a fit with a superposition of a central 
line and eight satellites (see text). Their resonance position has been calculated with the Breit-Rabi 



formula (Ref. 



26 



23) with a hyperfine constant A( 167 Er)=71(l) Oe. 



The temperature dependence of the linewidth obtained from the fits are displayed in 
Fig. [21 At low temperatures we observe a linear increase of the linewidth with temperature 
(Korringa relaxation). Towards higher temperatures the linewidth grows faster than linearly 
and above ~ 20 K the linewidth is of the order of the resonance fields and the EPR is not 
detectable any more. 

The non-linear temperature dependence of the EPR linewidth of non-S" ground state rare 
earth moments in metals can be conveniently described by the function 



AH = AH res + bT + 



(17) 



exp(A/kBT) — 1 

where AH res represents the residual linewidth, bT the Korringa relaxation of the local- 
ized moments to the conduction electrons and the supralinear third term summarizes the 
contribution from relaxation processes via excited CEF states located at an energy A above 
the CEF ground state (Orbach and Hirst processes).— 1 ^ Fitting our data with eq. (JI7|) we 
arrive at a Korringa term 
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FIG. 2: (Color online) EPR linewidth AH as a function of temperature for Lui-^Er^Gaa with 
x=0.005, 0.01 and 0.03. The supralinear increase towards higher temperatures can be fitted (dashed 
line) assuming an exponential increase according to eq. (|17p with an energy gap from the ground 
state to the first excited state of A = 35 K or 3.1 meV. 



The exponential broadening of the linewidth at higher temperatures can be ascribed to 
relaxation via an excited CEF state at an energy 



The g factor of the Er 3+ EPR in Lui_ x Er x Ga3 clearly confirms the Tj doublet as the 
CEF groundstate and the analysis of the temperature dependence of the linewidth sets the 
first excited state to ~3.1 meV, in very good agreement with the energy of the first excited 
quadruplet concluded by Murasik et ai. These findings strongly support the analysis by 
Murasik et al. for the CEF level scheme of Er 3+ in ErGa3. They also show that by dilution 
into the isotypic diamagnetic matrix, LuGa3, the single-ion CEF parameters remain close 
to those found in the concentrated compound in the paramagnetic state. 

EPR spectra taken at low temperatures on polycrystalline samples of Luo.95Tm .o5Ga3 
showed no indication of a resonance from Tm 3+ . 



b = 30(2)Oe/K 



(18) 



A = 3.1(5)meV 



(19) 
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B. Inelastic Neutron Spectroscopy 



In their early INS investigation Morin et al. observed broad inelastic modes centered at 
~2.6 meV, ~llmeV and ~17 meV from which they concluded a total splitting of the CEF 
levels of ~200K. 11 Guided by these results we studied polycrystalline samples of TmGa 3 , 
Luo.95Tnio.o5Ga3 and ErGa3 on ILL's time-of-flight spectrometer IN4 with neutrons of inci- 
dent energies corresponding to wavelengths of 3.3 A, 2.2 A and 1.1 A at temperatures of 20 
K and below. Our data for TmGa3 show modes at the same energies as found by Morin 
et al. but the mode at ~llmeV is somewhat better resolved and it becomes obvious that 
this mode consists of two modes at ~8.5meV and ~10.7meV. At the same energies, modes 
of comparable (normalized) intensity are also observed in the IN4 spectra of ErGa3 and of 
Luo.gsTmo.osGas. In ErGa3 Murasik et al. have seen inelastic scattering intensity emerging 
from CEF transitions essentially up to ~5meV.— This energy limitation for the CEF tran- 
sitions in ErGa3 and the fact that in the diamagnetically diluted sample Luo.gsTmo.osGas 
similar modes are observed, leads us to the conclusion that all modes appearing in TmGa3, 
ErGa3 and Luo.gsTmo.osGas above ~6meV cannot be attributed to CEF transitions. The 
subsequent inspection of the ^-dependence of these modes clearly identified them as phonon 
modes. Extending the energy range further, two more phonon modes at ~21.5meV and 
~35meV were found. 

The detailed investigation of the temperature dependence of the low energy regime, 
Ei <7.5meV, showed that the mode centered at ~2.6meV consists at least of two sub- 
modes and revealed also a mode at very low energy <0.4meV appearing as a shoulder on 
the energy-loss side of the elastic peak. Comparison with the energy-gain part of the spec- 
tra indicated this mode to originate from a CEF excitation from states close to the ground 
state. Summarizing these results, we concluded that, in contrast to the preceding analysis 
by Morin et al. the total CEF splitting of the J=6 manifold of Tm 3+ in TmGa3 does not 
exceed an energy range of ~3.5meV. This presumption was in agreement with the energy 
range estimated from the A± and Aq parameters obtained from the analysis of the ErGa3 
INS data (see above). Using the IN4 data a least square fitting procedure was carried out 
on the 5 and 10 K spectra (above the phase transitions). Although the 0.3meV excitation is 
not clearly resolved, the relative intensities of the peaks at 0.3meV and 2.5meV enabled us 
to obtain a first estimate of the LLW parameters which were in accord with this conclusion. 
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In order to better resolve the energy regime below ~4 meV, we performed measurements 
using lower incident energies on the cold source TOF spectrometer IN6. To further improve 
the resolution in the regime below ~1 meV, spectra of the diluted system Luo.9Tmo.1Ga3 
were taken in the cold neutron regime on ILL's triple-axis spectrometer IN12. Fig J3] displays 
a typical INS spectrum of TmGa3 measured on IN6 at 5K in the paramagnetic regime and 
of Luo.9Tmo.1Ga3 obtained with the triple-axis spectrometer IN12 at 2K. The spectrum 
of TmGa3 clearly shows a broad structured mode centered at 2.6meV with a FWHM of 
~0.8meV, larger than the experimental resolution. This broadening indicates that this 
mode consists of several overlapping lines, as will be discussed in detail below. A strong 
mode overlapping with the elastic line on the energy-loss side is now well resolved. It 
becomes more clearly visible as a mode separated from the elastic line in the IN12 spectra 
of Luo.9Tmo.1Ga3 at 2K. 

To fit the INS spectra we used the program Mcphase22- which provides the INS total 
scattering cross section (in barn/ion) for a transition from level n to m at a set temperature 
with the partition function Z and with LLW parameters xllw and W giving the CEF 
energies Ei according to eq. ( 1201) . 

V* - 4, (g) W-W) x 2 E ! < n]Am > |2 (20) 

An additional program was set up to fit the spectra using a wide range of xllw and W 
using a Lorentzian lineshape broadening of the modes. To fit the spectral range close to 
origin we also included an elastic peak. 

Using the data collected on IN6 for six temperatures 5, 10, 15, 20, 30 and 50 K and using 
the estimate for xllw and W as obtained from the IN4 spectra as starting values, we fitted 
the spectra using the intensities and energies listed in Table [H Good fits over the whole 
temperature range up to 50K were obtained for all temperatures with 

x LLW = -0.44(1) and W = -0.222(2)K. (21) 

The fits for all temperatures taken on IN6 are shown in fig HI Due to there being many 
transitions very close to each other (see tabled]) the modes at 2.37 and 2.41 meV were merged 
into one mode. As well as the peaks at 2.59 and 2.61 meV and also the levels at ±0.34 and 
at ±0.40 meV. The transition at 3.03 meV is increasingly intense at lower temperatures 
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FIG. 3: INS spectra (o) of (a) Luo.9Tmo.1Ga3 collected on the triple-axis spectrometer IN12 at 2K 
(the solid line is a guide to the eye) and of (b) TmGa3 collected on the TOF spectrometer IN6 at 
5K. 

(see table [I]) and was included into the fit at 5 and 10 K. Figure H] gives the comparison of 
experiment and calculated spectra. 

Further analysis of the transitions shows that at 30 and 50 K the agreement with the 
fit profiles and the measurement is always better than at lower temperatures. Especially 
at lower temperatures it became increasingly difficult to fit the energy range between the 
transition at 0.4 meV and 2.5 meV. One can see that at 20K the agreement with the data 
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TABLE I: Transition probabilities according to ([20]) for x LLW = -0.44 and W = -0.222 K of the 
crystal field states of Tm 3+ in TmGa3 



r —> r 


A E (meV) 


5K 


10K 


Urn ( ba] 
15K 


•n/ion) 
20K 


30K 


50K 


I2 -> i 5 


0.4 


7.14 


5.07 


4.30 


3.86 


3.37 


2.94 


ri - r 4 


0.34 


6.75 


5.37 


4.73 


4.33 


3.85 


3.41 


P (2) _^ r 


-0.4 


2.85 


3.20 


3.17 


3.07 


2.89 


2.68 


r 4 -> ri 


-0.34 


3.08 


3.63 


3.64 


3.56 


3.38 


3.15 


rf } - r 3 


2.41 


1.35 


1.52 


1.50 


1.46 


1.37 


1.27 


rf } - il 


2.63 


0.79 


0.88 


0.87 


0.85 


0.80 


0.74 


r 4 - r 3 


2.37 


0.44 


0.52 


0.52 


0.51 


0.48 


0.45 


J-4 -> 1 5 


2.59 


0.91 


1.07 


1.08 


1.05 


1.00 


0.93 


1 2 -> 1 5 


3.03 


1.34 


0.95 


0.8 


0.72 


0.63 


0.55 



and the fit is good for the transition energies and the energy regions around them, but for 
the part of the spectrum between 1 meV and 1.5 meV there is a gap between the fit and 
the data, this persists from 20 to 5 K. The reason for this discrepancy may be due to a 
quasielastic contribution which one would expect to decrease with increasing temperature. 
A quasielastic line was not included in the fits. Additionally, at these low energies, there 
could also be some broadening due to magnetic dispersion. 

The 50 K spectrum in comparison with the fit is displayed in more detail in fig. [5j The 
calculated energy positions of the modes according to Table [J are also indicated. The final 
CEF level scheme with the parameters given in eq. (122]) is shown in fig [6] with a comparison 
to the CEF level scheme proposed by Morin et al..— The most intense transitions (T n ^ T m ) 
(cf. Table [I]) are shown for energy gain and loss. 

In order to reduce magnetic and quadrupolar and exchange interaction between the Tm 
moments we studied the diluted Tm sample, Tmo.1Luo.9Ga3 . Heat capacity measurements 
on these sample showed no magnetic or quadrupolar ordering down to 1.8 K (see below). 
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FIG. 4: (Color online) Semi-log plot of the IN6 spectra taken at 5, 10, 15, 20, 30 and 50 K. Black 
solid line represents that fit as described in the text. The green solid lines are the individual 
Lorentzian profiles. 




-0.5 0.0 0.5 1.0 

energy transfer (meV) 

FIG. 5: (Color online) 50K spectra of TmGa3. Crystal Field transitions are shown for xllw 
-0.1 1 and W = -0.222 K for the low excitations. Insert: 2.5 meV peak shown to consist of 4 modes. 
The red line is the fit obtained by fixing the intensities and positions obtained from the program 
Mcphase and using a Lorentzian profile for each line. 
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FIG. 6: left: Level scheme proposed by Morin et al.M Right: Our level scheme, note that the 
ground state (T2) has been shifted down to -0.2 meV so to be clearly seen from the T\ which is 
only 0.009 meV away. Solid arrows represent energy gain transfer, dash dot arrows energy loss. 
Only the most intense transitions have been shown and for T n ^ T m . 

The spectrum at 2 K around the elastic peak measured on IN 12 on the diluted sample 
is shown in fig [7J The insert displays the full spectrum at 2 K. The transitions at E ~ 
±0.3 meV can be clearly resolved being well separated from the elastic peak. The mode at 
~2.6meV is also seen but it is as broad as for the concentrated TmGa3 sample. Additionally, 
there is a small feature at ~1 meV which is not seen in the other data on IN4 on the Tm 
concentrated compound TmGa3 and IN6. Data taken at 10 K also showed similar features 
as those at 2 K, including the peak at 1 meV. We tentatively ascribe this extra feature to a 
splitting of the Y\ third excited state. Such a splitting is probably caused by some slight 
local deviations from cubic symmetry due to the random substitution of the Lu sites by Tm 
atoms in the diluted system Tmo.1Luo.9Ga3. Such a splitting could also be the reason why 
the mode at ~2.6meV is as broad as in the concentrated sample. By inspecting the relative 
intensities and also the energy positions, we find slight deviations in xllw and W from the 
pure TmGa3 to the diluted case, which can be expected because of the slightly different 
lattice metrics of LuGa3 with respect to that of TmGa3. Such a slight reduction of the W 
parameter for single ion CEF splitting of Tm in LuGa3 and TmGa3 is also confirmed by the 
analysis of the specific heat data (see below). 
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FIG. 7: (Color online) Semi-log plot of the IN12 data for Tmo.1Luo.9Ga3 at 2 K. Open black circles 
are the data points, solid line represents the fit as described in the text. Insert: Full IN12 spectra 
from -0.6 meV to 3.5 meV. 

The best fit of the CEF spectra of Tnio.1Luo.9Ga3 is obtained with the parameter set 



^LLW 



-0.44 and W 



-0.18K. 



(22) 



C. Specific Heat and Magnetisation Measurements 

Our INS results indicate a CEF scheme for TmGa3 which is significantly different from 
that proposed by Morin et ai. A characteristic feature of our CEF scheme is a separation 
of the CEF levels into two groups about ~2 meV apart. The low energy group contains 4 
CEF states (r2 Ti, and with an overall splitting of less than ~5K. The splitting 
of the CEF levels in the group around 2 meV (r3 and T§) is even smaller and amounts 
to < 2K. The small overall excitation energies of the CEF level scheme of ~25 K allow 
us to observe their contributions, C m , to the total heat capacity at low temperatures and 
to clearly distinguish these magnetic contribution from phonon contributions which, at low 
temperatures, decrease rapidly according to C p h oc T 3 . In order to determine the mag- 
netic contribution to the heat capacity, we measured the heat capacity of samples with the 
composition Lui_ x Tm x .Ga3 (x=0.05 and 0.1) in the temperature range 1.8 - 40 K and of 
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TmGa3 in the range 2 - 100 K and compared their heat capacities with the heat capacity 
of LuGa3, the latter representing the phonon contribution. Samples with Tm diluted into 
the diamagnetic host LuGa3 allow to isolate - within limitations of cluster formation with 
concomitant intercluster Tm - Tm interaction - the single ion CEF behavior and to follow, 
by increasing the Tm content, the effects of exchange and quadrupolar interaction. 

Fig [8] shows the resulting difference heat capacities scaled by 1/0.05, 1/0.1 and 1 for 
Luo.gsTmo.osGas, Luo.9Tmo.1Ga3 and for TmGa3, respectively. In order to compensate for 
small mismatches of the high temperature data between Luo.95Tm .o5Ga 3 and LuGa3, the 
heat capacity of LuGa3 was adjusted by the factor 0.995 such that at high temperatures the 
resulting difference magnetic heat capacities fall off oc 1/T 2 . Also shown in Fig [S] are the 
magnetic heat capacities calculated from our CEF level scheme and for that given by Morin 
et al..— In contrast to the expected scenario according to the CEF scheme proposed by Morin 
et al. our data reveal essential magnetic contributions only up to about 30 K. They split into 
two characteristic anomalies: A broad Schottky-type contribution, with the maximum at 
~10K, followed by a sharp increase towards lowest measuring temperatures. These features 
are very well reproduced by the calculations using the CEF level scheme derived from our 
INS measurements. By doubling the Tm concentration i.e. going from Luo.gsTmo.osGas 
to Luo.9Tmo.1Ga3, the dip around ~ 5 K between both anomalies is somewhat smeared 
out and the low temperature increase becomes slightly more pronounced. TmGa3 shows 
a very sharp A-like anomaly at the AFQ and AFM ordering temperature near ~4.27(1)K 
and a broad Schottky-like anomaly centered at ~10K, very similar to those observed for 
the diluted systems but somewhat more smeared, with a long tail cx 1/T 2 clearly visible up 
to temperatures >35K. It appears that the magnetic entropy gained in the diluted systems 
below ~5K for TmGa3 merged into the very sharp anomaly near the long-range ordering 
point. 

Simulating the temperature dependant susceptibility and the magnetisation provides an- 
other way of testing the CEF parameters. Fig. [9] shows a plot of the inverse susceptibility 
in the range of 2 K to 200 K. The susceptibilties obtained from the CEF parameters as 
obtained by our inelastic neutron spectroscopy agree well with the experimental data above 
the transition temperatures. Differences of the simulated susceptibilities using the CEF pa- 
rameters xllw an d W given by Morin et al.— to the experimental data appear below ~50 
K. The good agreement of the fit lends further support to our CEF scheme. 
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FIG. 8: (Color online) Magnetic contributions (per mole Tm atoms) for Luo.gsTmo.osGas (black 
□), Luo.9Tmo.1Ga3 (red o) and TmGa3 (blue A) obtained as described in detail in the text. The 
(black) solid line represents the magnetic heat capacity calculated for Tm 3+ in a cubic crystal 
electric field parametrized by the Lea Leask Wolf parameters xllw=-0-44 and W=-0.18K; the 



dotted (red) line t 
Morin et al. (Ref. 



le magnetic heat capacity calculated for ccllw=-0-32 and VF=1.03K given by 
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and the (green) dashed curve represents the heat capacity of LuGa3 taken 



as reference for the lattice heat capacity. 

The magnetisation of TmGa3 in the paramagnetic phase with magnetic field applied 
along some main symmetry directions was also calculated in order to check for consistency 
with our CEF model. The results are shown in fig [TU] in comparison with data collected at 
10K by Morin et al. along [001] and [110] directions.— We also display the magnetisations 
calculated using the CEF scheme proposed earlier by Morin et al. which leads to significant 
discrepancies. In contrast, our CEF scheme provides an almost perfect description of the 
magnetism of TmGa3 in the paramagnetic phase. 



D. Effect of a tetragonal distortion on the CEF level scheme 

The analysis of the heat capacity and the magnetisation of TmGa 3 and the samples 
Lui-^Tm^Gas provide very strong evidence that the parameter set xllw= -0.44 and W= 
-0.22 K obtained from the INS measurements describe the CEF scheme of Tm 3+ in TmGa 3 
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FIG. 9: Temperature dependence of the inverse Susceptibility of TmGa3. Open circles represent the 
data. Solid black line the calculation using the new crystal field parameters. Dash line calculation 
with the parameters of Morin et al.. 31 

correctly. We now discuss possible effects of the quadrupolar ordering on the CEF level 
scheme. Quadrupolar ordering lowers the cubic symmetry and one has to consider additional 
terms B^O® and B\0\ in the CEF Hamiltonian (see eq. If one continuously turns on 

the S° term a splitting of the levels as shown in fig [TT] results. The low energy states become 
mixed already for small values of B\ being of the order of B\ = -0.3 K as reported Gubbens 
et al. based on Mossbauer spectroscopy data.— 

V. DISCUSSION 

Our new set of crystal field parameters as obtained from inelastic neutron scattering give 
thermodynamic properties which agree very well with the experimental data. We have shown 
that the CEF scheme proposed by Morin et al. disagrees with thermodynamic properties 
and should be discarded. As the crystal field energy levels in TmGa3 are very close to each 
other, and are all below 4 meV, determining such a level scheme is difficult. To ensure 
that we have the correct parameter set, we have performed a number of tests. First, using 
the results of Murasik for ErGa3, by EPR measurement on doped Er samples, we have 
confirmed that the level scheme by Murasik et al. is correct. We have then scaled these 
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FIG. 10: Magnetisation at 10 K along Top: [110] Bottom: [001]. Data taken from Morin et al. 
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values to give us an approximation of where in the parameter space the correct values of 
xllw and W for TmGa3 could be. As neutron spectroscopy allows the best determination of 
such parameters, we have fitted numerous spectra and determined values that best suit all 
spectra. These values have been taken to calculate thermodynamic properties, such as the 
Schottky anomalies in the specific heat and the temperature dependant susceptibility as well 
as the magnetisation in several lattice directions. The agreement with the values obtained 
by inelastic neutron spectroscopy and the thermodynamic properties assures us that we have 
obtained the correct parameter set. Due to the small energy range of the whole CEF level 
scheme AFQ ordering in TmGa3 must involve an essential mixing of the CEF levels. Our 
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FIG. 11: Effect of a B5> term in the Hamiltonian eq. ([!]) on the CEF level scheme of TmGa3. 

results can be used for further measurements in the ordered phase in order to understand 
the complex and interesting competition between orbital and magnetic ordering in TmGa3. 
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